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Abstract 

One-loop off-shell contributions to the quark-gluon vertex are calculated, in 
an arbitrary covariant gauge and in arbitrary space-time dimension, including 
quark-mass effects. It is shown how one can get results for all on-shell limits 
of interest directly from the off-shell expressions. In order to demonstrate that 
the Ward-Slavnov-Taylor identity for the quark-gluon vertex is satisfied, we have 
also calculated the corresponding one-loop contribution involving the quark-quark- 
ghost-ghost vertex. 
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1 Introduction 



Tests of perturbative Quantum Chromodynamics 0,0] are steadily reaching new levels 
of precision , necessitating information on higher- loop results for a wide range of pro- 
cesses. Higher-order QCD effects are also frequently required for background estimates 
in searches for signs of new physics. 

In spite of its fundamental role, the quark-gluon vertex has not been explored in 
detail, even at one loop. It is of course very much related to the electron-photon vertex 
of Quantum Electrodynamics, a basic and nontrivial aspect of which is the anomalous 
magnetic moment, or g — 2, which provides a powerful test of the whole concept of 
Quantum Field Theories. The quark-gluon vertex differs from the electron-photon vertex 
already at one loop, by the contributions of an additional Feynman diagram, involving 
the three-gluon vertex. In fact, apart from introducing additional colour structure, this 
non-Abelian diagram introduces at the one-loop level a kinematical structure which is 
absent in the QED vertex. 

For special kinematical configurations, and special gauges, several results are avail- 
able. Already around 1980, the symmetric off-shell case was considered by Pascual and 
Tarrach (PT) in [^] (see also |^) in an arbitrary covariant gauge for massless quarksQ. 
The emphasis was on comparing the MS and Weinberg's renormalization schemes. The 
symmetric off-shell case was also considered in ||^ by Dung, Tarasov and Phuoc (DTP), 
for massive quarks, restricted to the scalar function multiplying 7^. For massless quarks, 
some on-shell results are available, mainly in the Feynman gauge, presented by Nowak, 
Praszalowicz and Slomihski (NPS) §]. 

The situation when one (gluon or quark) momentum vanishes has been studied in 
more detail. Technically, in this case all three-point functions effectively reduce to two- 
point integrals. For massless quarks, some one-loop results in arbitrary covariant gauge 
have been obtained by Braaten and Leveille ||^. In the Feynman gauge, also two- loop 
corrections have been presented in 0. Moreover, in a recent paper by Chetyrkin and 
Retey ||10[, renormalized expressions for three-loop-order QCD vertices have been ob- 
tained for such zero-momentum configurations, in an arbitrary covariant gauge. 

The QED contribution, proportional to the "Abelian" contribution to the quark- 
gluon vertex, has been studied more systematically. An early paper by Ball and Chiu 
(BC) presented a systematic kinematical decomposition of the vertex, and gave off- 
shell results for the one-loop QED vertex in Feynman gauge. Their work was extended 
to arbitrary covariant gauge by Kizilersii, Reenders and Pennington (KRP) [Q. All 
above-mentioned papers deal with the (dimensionally-regulated |TB]) four-dimensional 



case. Results for the three-dimensional QED contribution are also available (for massless 
fermions), due to Bashir, Kizilersii and Pennington (BKP) [|T^. A summary of all these 
one- loop results is given in Table 1. In addition to this table, we note that another 
special gauge which has been investigated is the Fried- Yennie gauge |]I5 |. 

Among non-covariant gauges, we would like to mention the Coulomb gauge. In 
some sense, it is more "physical", but technically rather challenging |T^. A rather 



^Numerical analysis of two-loop QCD vertices in this limit is given in a recent paper 
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Table 1: Kinematics and gauges considered in other one- loop studies. None of these 
results is valid for arbitrary dimension n. 



different approach to QCD vertex functions is provided by lattice calculations [T^]- The 
quark-gluon vertex functions may also serve as a basis for modeling the photon-nucleon 



vertices and the quark-Reggeon vertex 119 



From Table 1, one can see that, even if we consider the results in (or around) four 
dimensions, there are still several "white spots". The aim of the present paper is to 
cover all such remaining spots. Moreover, we present results which are valid for an 
arbitrary value of the space-time dimension. Apart from the quark-gluon vertex itself, 
we also consider the related two-point functions, and the quark-quark-ghost-ghost vertex 
function, in order to be able to check that the obtained quark-gluon vertex function obeys 
the Ward-Slavnov-Taylor (WST) identity [g^]. 

At the one-loop level, the simple and well-known Dirac-matrix structure of the lowest- 
order quark-gluon vertex gets modified. In the general case, twelve structures are needed 



to decompose it |TT|. Thus, twelve scalar functions multiplying these tensor structures 
are to be calculated. These scalar functions depend on the gauge parameter, the space- 
time dimension, quark mass(es), and the kinematical invariants {pf, p\^ p^. Four of 
them (the "longitudinal" ones) are involved in the WST identity, whereas the remaining 
eight are unconstrained. 

There are several reasons why the one-loop results calculated in arbitrary gauge and 
dimension n are of special interest: 

(i) knowing the results in arbitrary gauge, one can explicitly keep track of gauge invari- 
ance for physical quantities; 

(ii) if one is interested in the two-loop calculation of the quark-gluon coupling, one should 
know one-loop contributions in more detail; 

(iii) results in arbitrary dimension make it possible to consider all on-shell limits directly 
from these expressions (see Section 4), this is impossible if one only has the results valid 
around four dimensions; 

(iv) QCD is also a theory of interest in three and two dimensions (see, e.g., pl|); 



(v) as we shall see, the results for arbitrary dimension are not much more cumbersome 
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than those considered around four dimensions (in some respects, they are even more 
transparent and instructive). 

The paper is organized as follows. In Section 2, we introduce the notation for the two- 
and three-point functions to be considered, and discuss their decomposition in terms of 
scalar functions as well as the corresponding Ward-Slavnov-Taylor identity. In Section 3, 
we present the most general off-shell results for the quark-gluon vertex. Section 4 contains 
the corresponding expressions for special limits of interest. In Section 5, we conclude 
with a summary and a discussion of the results. Then, we have several appendices where 
some further results and technical details are presented, such as the formulae used to 
decompose the quark-gluon vertex (Appendix A), relevant results for the scalar integrals 
(Appendices B and C), results for the one-loop contribution involved in checking the 
WST identity (Appendix D), and general results for the transverse part of the quark- 
gluon vertex (Appendix E). 

2 Preliminaries 

We shall here establish some notation, and discuss the functions involved in the Ward- 
Slavnov-Taylor identity for the quark-gluon vertex. 

2.1 Notation 

A graphical representation of the quark-gluon vertex is given in Fig. 1.0 The momentum 




Figure 1: Kinematics of the quark-gluon vertex. 

of the outgoing quark is denoted by pi, p2 is the momentum of the incoming quark, 
whereas and /i are the momentum and the Lorentz index of the gluon, respectively. 
All momenta are ingoing, pi + P2 + P3 = 0. The lowest-order quark-gluon vertex is 

9inM,^, (2.1) 

■''To produce the figures, the AXODRAW package |2|] was used. 
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where are colour matrices corresponding to the fundamental representation of the 
gauge group. As a rule, it will be implied that the SU(A^) group is considered, with 
being the number of colours (we can put = 3 in the end). 

When one calculates radiative corrections to the quark-gluon vertex, other Dirac 
matrix structures arise, in addition to 7^, Eq. ( |2.1| ). The total number of independent 
structures is 12 (see, e.g., in and also in Appendix A of this paper). Extracting the 

over-all colour structure, we can present the one-particle irreducible quark-gluon vertex 
as 

r^(pi,P2,P3) = ^?T'^r^(pi,P2,P3) , (2.2) 

where matrix notation in both colour and Dirac matrices is understood. 

At the one-loop level, we have two contributions to the quark-gluon vertex which 
are shown in Fig. 2. Their colour factors are proportional to {Cp — \C a) and Ca-, 
respectively, where Cp and Ca denote eigenvalues of the quadratic Casimir operator in 
the fundamental and adjoint representations, respectively. For the SU(A^) gauge group, 

N"^ - 1 

Ca = N, Cp = . (2.3) 




-Pi - Q 




\P2-(1 



Figure 2: The two one- loop diagrams. 



The first, "Abelian" contribution is completely similar to the one- loop correction to 
the fermion-photon vertex in QED. The difference is only in the over-all factor. Formally, 
we can get the one-loop QED vertex from the considered QCD vertex by putting Ca = 0. 
The second diagram in Fig. 2 is essentially non-Abelian and appears due to the self 
interaction of gluons. 

If quarks are massive, it is clear that diagrams a and b in Fig. 2 involve "triangle" 
integrals with two and one massive lines, respectively: 

JM, ^2, z^a) = J [(p,_,)2]^i [{p, + q)T [q'-mT' ^^'^^ 
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where n = A — 2e is the space-time dimension in the framework of dimensional regu- 



larization |T^. Understanding the subscript of J as the number of massive propagators 



(cf. in |23]), we can extend this notation to the massless integrals J(z/i, z/2, z/3) considered 
in H, via 

J{ui, U2, 1^3) ^ M^i^ ^2,1^3) ■ (2.6) 

Integrals with Lorentz indices can be reduced to the scalar ones using the standard 
techniques p5| (see also in p6|,^). Using the integration-by-parts technique p8| (see 
also in PS), all integrals with higher integer powers of propagators can be algebraically 



reduced to integrals with the powers equal to one or zero (for details, see Appendix B) 



As in , we shall extract from the expressions for one-loop integrals a factor 



V = 



1) 



T{n - 3) 



r(3 



T{l-2e) 



(2.7) 

^1,2) 



A natural extension of the notation used in ||2^ is to introduce the functions (pi {i 
such that 

J,(l, 1, 1) = i 7r"/2 ^ ^.{pIpIpI m). (2.8) 

In this sense, the function Eq. (2.14) of ||2^, would correspond to yjo, which also 
represents the massless limit of {i = 1,2). Moreover, we can reserve the notations 
J3 and (^3 for the triangle integral with all three massive lines, which occurs in the 
three-gluon vertex (the massive quark loop contribution). 
Then, for the two-point integrals we introduce functions 



(2.9) 



where pi [l = 1,2,3) is the external momentum of the two-point function, whereas the 
subscript {i = 0,1,2) shows how many of the two internal propagators are massive. In 
this way, kq (coinciding with the k defined in Eq. (2.15) of |2^) corresponds to the 
two-point function with massless lines. 



Ji(l,l,0) = Jo(l,l,0) = i7r"/2r/«:o,3 



(2.10) 



and analogously for Jo(0, 1, 1) and 70(1,0, 1), with ko,i and kq^2, respectively. Then, ki 
corresponds to the two-point function with one massive and one massless line. 



Ji(0,l,l) = J2(0,l,l) = i7r'^/'r//ci,i, 
Ji(l,0,l) = J2(l,0,l)=i7r"/2r//€i,2. 
Finally, K2 corresponds to the two-point function with two massive lines, 

J2(l,l,0) = J3(l,l,0) = i7r"/2^/€2,3, 



(2.11) 
(2.12) 

(2.13) 



and similarly for 73(1, 0,1) and J3(0, 1,1) (which would involve fi;2,2 and ^2,1, respec- 



tively). The massless two-point functions introduced in Eq. (2.15) of p4[ can be identified 
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The new feature, as compared to the massless case, is the appearance of the "tadpole" 
integral 

Ji(0, 0, 1) = J2(l, 0, 0) = J2(0, 1,0) = m^l'' ^^^(^')'~' = ivr"/' r^m^K, (2.14) 



with 



Let us recall that the massless tadpoles vanish in the framework of dimensional regular- 



ization 



Ji(l, 0, 0) = Ji(0, 1, 0) = J2(0, 0, 1) = 0. (2.16) 

We shall also introduce some notation to keep track of the various orders in the 
perturbative expansion. For a quantity X (e.g. any of the scalar functions contributing 
to the propagators or the vertices), we shall denote the zero-loop-order contribution as 
X(°), and the one-loop-order contribution as X'^^^ so that the perturbative expansion 
looks like 

X = X(°)+X« + ... (2.17) 

2.2 Two-point functions 

The lowest-order gluon propagator is 

- W'^^'^^ ^ (w. - e ^) , (2.18) 

where ^ is the gauge parameter corresponding to a general covariant gauge, defined 
such that ^ = is the Feynman gauge. Here and henceforth, a causal prescription is 
understood, 1/p^ ^ l/(p^ + iO). For the present purposes, loop corrections to the gluon 
propagator (see, e.g., Eqs. (2.7) and (C.l) of |^) are not required. 

We shall denote the quark propagator as S{j)). The two scalar functions a(p^) and 
/3(p^) in the inverse quark propagator are defined via 

\S-\v) ^ ^ fi{j?)l . (2.19) 

where j) = p^7^, whereas I is the unit matrix in the space of Dirac matrices. At the 
lowest order, a^^'^ = 1 and P^^^ = —m. For the next-to-leading order, one needs to 
calculate the one-loop diagram shown in Fig. 3, which yields (for n near 4, see [pO[] ) 

a^'\p') = j£^^{n-2)il-0[ip' + rn')K,{p';m)-m'K{m'^^ , (2.20) 
P^'\p') = m{n-i)K,{p'-m) . (2.21) 
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P p-q P 



Figure 3: Quark self-energy diagram. 



The ghost propagator is 

^aia2 ^p2^ ^ • ^aiaa ^(pj ^ (2.22) 



p2 



The lowest-order result is G^'^^ = 1, whereas the one-loop contribution reads 

G^'\p') = % [2 + (n - m] n,{p\ (2.23) 

Note that in the Fried- Yennie gauge (see also in ^3), ^ = —2, G'*^"'^)(p^) is finite as 
n — > 4. Moreover, if one chooses ^ = —2/ (n — 3) as the n-dimensional generalization of 
this gauge ||33|,0, then the r.h.s. of Eq. ( |2.23|) vanishes. 



2.3 Ward— Slavnov— Taylor identity 

The WST identity for the quark-gluon vertex r^(pi,p2,P3) reads (see, e.g., in |]^,|3^ 

P'z ^M,P2,P3) = G{pI) [S-\-pi) Hip,,p2,Ps)~Hip2,Pi,P3) S-\p2)] , (2.24) 

where G{p'^) (see Eq. ( |2.22D ) is the scalar function associated with the ghost propagator. 

The function H (and the "conjugated" function H) involves the complete four-point 
quark-quark-ghost-ghost vertex, as shown in Fig. 4. To get the H function, we need to 
"join" the out-quark and out-ghost lines in a non-standard vertex (denoted by a cross) 
and integrate over the resulting loop momentum. It should be noted that the complete 
quark-quark-ghost-ghost vertex involved in the WST identity can be decomposed into a 
connected and an unconnected piece, as shown in Fig. 5. Moreover, the connected part 
can be further split in terms of the proper (one-particle irreducible) vertices, see Fig. 6. 

We note that the first diagram on the r.h.s. of the equation shown in Fig. 6 (the 
diagram involving the proper four-point function) does not not have a zero-loop (tree) 
contribution. Its perturbative expansion starts from the one-loop boxes, corresponding 
to the exchange by two gluons. Since the H function involves an extra loop integration 
(see Fig. 4), this proper four-point function does not contribute to the one- loop-order H 
function, H^^\ which is shown in Fig. 7. 
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k 



Figure 4: Graphical representation of the H function. 




unconnected part connected part 



Figure 5: Unconnected and connected parts of the quark-quark-ghost-ghost amphtude. 

The H function can be decomposed in terms of scalar functions ( "form factors" ) as 

(2.25) 




Figure 6: Connected part of the quark-quark-ghost-ghost amplitude in terms of proper 
(dark) vertices. 
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Figure 7: The one-loop order function H'^^\ 



with 

(^^lu = I ht^^u - 7i^7m) • (2.26) 
The "conjugated" function H can be written in terms of the same scalar functions, 

H{p2,Pi,P3) = Xo{phpl,pl)I - X2{pI,pI,pI)^i - Xi{phpl,pl)h + X3{pl,pl,pl)cr,iuPiP2- 

(2.27) 

At the lowest order, Xo'^ = 1 ^-^id xf'^ = (z = 1,2,3). The one-loop results for the Xi 



functions (valid for arbitrary values of n and ^) are presented in Appendix D. 

At the one-loop level, it is convenient to "split" the WST identity into two separate 
identities, corresponding to the contributions of the two diagrams shown in Fig. 2. To 
do this, we need to rewrite the one-loop contribution to the r.h.s. of ( p.24| ) in terms of 
colour coefficients {Cp — ^Ca) and Ca, in analogy with the two contributions to the 
l.h.s. On the r.h.s., all one- loop contributions are proportional to Ca, except for the 
quark self energies, which contain Cp. Therefore, all we need to do is to represent this 
Cp as {Cp — \Ca) + ^Ca- In this way, we get two separate WST identities for the 
contributions of diagrams a and b in Fig. 2, 



P^r(^"npi,P2,P3) 
p!;rf)(pi,p2,P3) 



ICa 



Cp' 



s 



-1/ 



-Pl] 



S-\-p,) H^'\p,,P2,P;)-H''\p,,p,,p,) S~\p 



S^\P2) 



(1) 



(2.28) 

(0) 



+\CaC~p' \S'\-pi 



+2G^'\pI) 



S' 



-Pl] 



S-\P2) 
S-\P2] 



(1) 



(0) 



(2.29) 



where, following the convention of ( 2.17) ), the superscripts "(0)" and "(1)" correspond 
to the zero-loop and one-loop contributions, respectively. 

The first identity, Eq. ( p.28|) , has, up to a factor, the same form as the Abelian (QED) 
identity, also known as the Ward-Fradkin-Takahashi identity |^6[|. The second identity, 
Eq. ( p.29| ), is the non- Abelian one. 
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2.4 Decomposition of the quark-gluon vertex 



Keeping in mind the WST identity ( p.24| ), it is useful to split the quark-gluon vertex 
into a longitudinal part and a transverse part, 



where 



P3'^^P{PUP2,P3) = 



(2.30) 

(2.31) 



and, therefore, F^J^ does not contribute to the l.h.s. of Eq. ( |2.24|) . 

In general, we shall just extend the decomposition of the QED vertex suggested 
in |11| (see also in Ref. [|12|) to the QCD case. The longitudinal part of the vertex can 



be represented as 



i=l 



with 



(2.32) 



(i^i -h)ipi -P2)/., 

{Pl - P2)m, 

(^f^iy {pi - P2y- 



(2.33) 



Using Eq. ( p.24D , the functions Aj can be related to the functions a, P, G and Xi- For 
instance, in the "Abelian" case (i.e., when we consider only the first, QED-like diagram 
in Fig. 2) the functions Aj {i = 1, 2, 3, 4) would be equal to (up to a colour factor) 



a{pl) + a{p. 



P1-P2 



and 0, 



(2.34) 



respectively (see in |lTT|, p!^j )p|. 

The transverse part of the vertex, which does not contribute to the WST identity 
( p.24| ), can be presented as |Tl[ 



^''P{Pi,P2,P3) =J2'^i(Pl^P2^P3) Ti,^,{pi,P2) 



(2.35) 



2 = 1 



^Since in the QED case the longitudinal functions have such simple representations (2.34) in terms of 
a{pi 2) and l3{pi 2)7 there was no need in Refs. to introduce a special notation for these functions. 

In the presence of the non- Abelian contribution, the situation becomes more complicated. This is why 



we have introduced longitudinal functions Xi in Eq. (2.32) 
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where the transverse tensors are the following: 

T2,^ = - [Pli,{P2P3) -P2^{PlP3)] (A -h) 



3, 



Ta,^ = [Pl,i{P2P3) -P2,,{PlP3)]cruXPlP2 

Te^f, = if,{pl -pI) + {pi- P2)j3 



Tj,^ = - pI) [7m(J^1 - h) - (Pl - P2),.] - (Pl - P2),,C^,yXPlP2, 

Ts,^ = -lf^cruxplP2 + Pl^lh - P2iii\- (2.36) 

The connection of A's and r's with the naive decomposition basis is discussed in Ap- 
pendix A. 

Applying charge conjugation to the quark-gluon vertex, i.e., interchanging quark and 
anti-quark, the following relation is obtained (see, e.g., in Ref. [IT 



c r^(pi,P2,P3) = -rJ(p2,Pi,P3). (2.37) 

Interchanging the quark momenta (pi ^ P2) and using the fact that 

C 7m = -7J, (2.38) 
one finds that all and are odd, except for L4^^ and Tg^^: 

Li,f,{pi,P2) = -Ll^{p2,Pi), i = 1,2,3, 

Li,^,{pi,P2) = Ll^{p2,Pi), 

Ti,^iPuP2) = -T^^ip2,Pi), i = 1,2,3,4,5,7,8, 

n,^{puP2) =T;[^,{P2,Pi). (2.39) 
To satisfy Eq. ( p. 37 ), all A's and r's must be symmetric under the interchange of pi and 



pI, except A4 and tq, which are odd: 

>^i{pl,pl,pl) = MpI,pI,pI), i = 1,2,3, 
^4(^1,^2,^3) = -^iipliPiiPl)^ 
'riipl,pl,pl) = nipl,pi,pl), i = 1,2,3,4,5,7,8, 
reiplplvl) = -reiplplpD- (2-40) 

An important corollary of these relations is that in the case pf = P2 = p^ the A4 and tq 
functions must vanish, 

K{p\p\pI) = 0, r^{p\p\pl) = . (2.41) 
12 



Furthermore, in Ref. a modification of tlie basis ( |2.35| )-( ^.36| ) lias been proposed, 
whicfi lias an advantage in dealing with kinematical singularities^. Namely, the transverse 
part is represented as 



^'■P{PI,P2,P3) =Y.'^iiPl^P2^Pl) Ti,^iPl^P2), 



(2.42) 



i=l 



where 



4, /I 2 2 

Pi- Pi 



T- = T- 



2T4,^ - = pilpif, - - 7m(J^i - 1^2)] - 2(pi + P2)^^(ru\ P1P2 



(^7^4), 



(2.43) 



with 



^4 = 4(^2 - Pl)^4, = r7 + 2P3r4. 



(2.44) 

Moreover, as we shall see below, in the on-shell limit the following modifications of A2 
and A3 turn out to be useful: 



A2 = A2 + 1^3 T-2, 



A3 = A3 - \pIti . 



(2.45) 



3 OfF-shell results 



Before presenting results for the Aj and Xj functions, let us introduce the following nota- 
tion for the Gram determinants occurring in the denominators: 



/C 
Ml 
M2 



2 2 
P1P2 



{pmf , 

/ 2 2\/ 2 2\ I 2 2 

[p^ — m )[P2 — m j + m P3 , 

{pj - w?){pI - m^)pl + rn?{pl - plf 



(3.1) 
(3.2) 
(3.3) 



In fact, /C is a symmetric function oip\, p\ and p\. It can be rewritten as — ;jA(p^,p|,p|), 
where X{x,y,z) is the Kallen function (cf. Eq. (3.2) of p^)- Note that M.i can also be 
represented as 



Ml 



2 2 
P1P2 



2{pip2)m^ + m'' = /C + {P1P2) + m 



To distinguish between the contributions of the two one-loop diagrams in Fig. 2, we 
shall use the letters a and h: 



A 



(1) 



A^) + A 



(lb) 



-(1) 



(la) 



(16) 



(3.4) 



For the calculation we used the algebraic programming system REDUCE Further 
technical details can be found in Appendices A and B. 

^In Ref. |l2|, the notation (Ji and Si was used for what we call Ti and Ti. 
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3.1 Results for the longitudinal part of the vertex 

The general results for the longitudinal functions of the vertex are reasonably compact, 
even in a general covariant gauge and arbitrary dimension, and are given below for the 
two diagrams. 

3.1.1 Diagram a 



g^T] [Cf-^Ca) (n-2)(l-0 



(47r)"/2 



X 



x^''\pI,pI,pI) = 



pKpI + rri^)!^!,! + p\{pI + m^)«i,2 - (p? + pDrri^i^l , (3.5) 
g^v{CF-lCA) (n-2)(l-0 



(47r)"/2 ^pIpHpI-pI) 



X 



>^t\plplpl) = 



p^iPi + m )ki,i - Pi(p2 + m )/«i,2 + (Pi - P2)m k 
g^T] (Cf - ^Ca) {n - 



(47r 



m/2 



„2_„2 ('«1,1-«1,2) , 
Pi P2 



a^)(pLp^,pD-o. 

Note that diagram a does not contribute to A4^^ . 



, (3.6) 

(3.7) 
(3.8) 



3.1.2 Diagram b 



x^i'\pI,pI,pI) 



x?\pI,pI,pI) 



9'vCa 1 



(47r)'*/2 16IC 

+[2 + (n- 3)^](pi - plf [{piP2)<Pi + mVi + «o,3 
+ [2+ (n - 3)C](pi -P2) [P2«;i,2 - P?/ti,i + (PiP2)(fi;i,i - fi;i,2) 

-(n-2)(2-C)/Cm2 
1 



, /«1,2 P?+P2. 



Pi 



pi 



Pipi 



(3.9) 



[2 + (n- 3)^1^3(^1-^2) (PiP2)<^i + m ipi + Ko,3 



(47r)"/2 16/C(pf-pi) 

+ [2 + (n - 3)^]P3 P2«^i,2 - pj^l,! + (PlP2)(«l,l - «;i,2) 



+(n-2)(2-0/C 



p? + pI + — po 2~ 

1 5 — Ki 2 H ^ ^ m K, 



Pi 



Pi 



Pipi 



,(3.10) 
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(47r)"/2 87Wi 



+ {n - 3)^ -v\- + {pI - pI - "^^)«i,2 + Ps'^o.j 



Pi 



+ (n-2)e 



P1-P2 

2{pip2) + 



^1,2 
P2 , 



pIpI 



+ [4(n 



<A^i^2 f- 

P1-P2 



(3.11) 



^^^C'a ^m{pl - pI) 
(47r)"/2 16/CA^i 



n-3)pi[{piP2)+m^f(^i + (n - 3)A:(/ti,i + /ti,2) 



+ICpl(pi + in-3)[{pip2)+m'^] [(piP3)«;i,l + (P2P3)/«l,2+P3^0,3 

,2 



+ (n-2)/C 



2{piP2) + m^ 2~ , -^1 



P1P2 



Pi - P2 V Pi 



pI 



(3.12) 



We have checked that these expressions, together with the results for the two-point 
functions and for the Xi functions (given in Appendix D), satisfy the WST identities 
( p.28| )- (|2.29| ) at the one-loop level, for arbitrary values of n and ^. 

Results for the Aj functions in the Feynman gauge (^ = 0) can be easily obtained 
from the expressions presented above. We just note that 



y^'i^'iplplpl) 
>^T\pIpIpI) 



5=0 



5=0 



q'^vCa (n — 1)"^ 
(47r)"/2 2(pf-pi) 

0. 



Ki 1 - Ki 2) 



(3.13) 
(3.14) 



3.2 Results for the transverse part of the vertex 



In an arbitrary covariant gauge, the results for most of the functions are rather cum- 
bersome. The expressions can be made more compact by introducing certain linear 
combinations of the Kjj functions, namely 



(PlP3)/«l,l + (P2P3)tl,2 + P3^i,3, + /^1,2 " 2Ki,3, 



(3.15) 



where we should take, for the subscript of Kj^3, i = 2 for diagram a and i = for diagram 
b. Note that the basis ( 3.15| ) is over-complete: we have got five combinations of the four 
k's. This degree of freedom has allowed us to present the results in a more compact 
form. 

We have collected results in an arbitrary covariant gauge in Appendix E. Here we 
present results for the Tj functions in the Feynman gauge. 
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3.2.1 Transverse functions in Feynman gauge: diagram a 



rt^iplplpl) 



^=0 



^4^''\pI,pI,pI) 

rt\pl,plpl) 
rt\plplpl) 



^=0 



^=0 



(47r)"/2 W 



- {pm) ^2 + 1^1,1 + 1^1,2 - 2«;2,3 



(3.16) 



g^r](CF- lCA) 1 r , „ 2, 2(n-l) 
2(n - 1) 



'pl<f2 



(]?1]?2) - m (piP3)«l,l + (P2P3)«1,2 + P3K-2,: 



-(n - 2) — (ki^i + Ki^s - 2k) - (n - 4) + ^1,2) 



+ (n-2) 



P1P2 

m'^{piP2){pi +pI) 
P1P2 



- 2w? - (pi -]?2)^ 



Pi - pi 



(3.17) 



{Cf - \Ca) 1 i2in-ll . 

16^\~^^^^^^ i^'^'' 



<P2 



+ ^^"'^ ^\ pi -P2f [{pm) - m^] [(piP3)/«l,l + (P2P3)«1,2 +P3/«2,3 



-8(n - 2)(]9^ -m^)(p2 - m^)^2 - 2(pi -P2)^ 4(n - l)m^ -^3 
+4(n - 2) [(P1P2) - (ki,! + «;i,2 - 2^2,3) 
+ (n - 4)(pi - p2)^ (/ti,i + ^1,2) 
m2(piP2) 



9^2 



+(n - 2) 

+(n - 2)m^ 
0, 



1-— (Pl - P2) - «1,2) 

(PlP2)(P? +P2) 



P1P2 



- 2 



(«i,i + «i,2 - 2k) 



(n - 4)m</72, 



^2^ (C^ - iC^) 

(47r)n/2 

(^i - I^-^) P^-pU o.. 2 2, 2(n-l) 

Wr/- lei^P"^ — ^ 



(3.18) 
(3.19) 

(3.20) 



P3 



(P1P2) - rn" 



<P2 



16 



2{n - 1) 
+ {n-2) 



m^{PiP2){pl + pI) 



PipI 



- 2m^ - {pi - p2f 



Pl-Pl 



r?^\plplpl) 
rt\plplpl) 



-{n - 2) — + Ki,2 - 2k) - - 4) (ki^i + Ki^a) 

P1P2 

0, 



(3.21) 
(3.22) 



+ (PlP3)«l,l + (P2P3)«1,2 + pIk-2,3]- (3.23) 

The function r4^"'' becomes non-zero when ^ 7^ (see Appendix E). In fact, the results 
for T^^""^ and Tg^^^ in an arbitrary covariant gauge can be obtained from Eqs. ( ^.16| ) and 
( p.23| ) by changing the overall factors n and (6 — n) into (n — ^) and [6 — n + {n — 4)^], 
respectively. 



3.2.2 Transverse functions in Feynman gauge: diagram b 



rl''\plplpl) 



5=0 



rt\plplpl) 



€=0 



g'^rjCA {n — l)m 



(47r)"/2 4/C 

-P1-P2 2 2- 

g^riCA 1 [2(r2-l) 



2[(PlP2) + - Ki^i - Ki^2 + 2^0,3 



(3.24) 



5=0 (47r)"/2 16/cj' /C 
2(n - 1) 



-Ps 



{P1P2) +m'^ V^i + 4(n - 3) {P1P2) + 



^1 



+2P3<^1 + 



(P1P2) + (PiPs) 1^1,1 + {P2P3) 1^1,2 + pIi^o,3 



-2(n - 3) (ki,i + Ki,2 - 2^0,3) + (n - 4) (/ti,i + Ki,2) 

-(^ - 2) .0 + /«i,2 - 2/t) - (n - 4)(pi - P2) —2 

P1P2 Pi ~ P2 



{PiP2){pI+pI) 



pIpI 



- 2 



pI-pI 



(3.25) 



P3{Pi-P2) {PiP2)+m 



(47r)"/2 32/C [ /C 
2(n - 1) 



(Pi - P2) (pm) + m {pm)ni^i + (P2P3)'«1,2 + P^i^o- 



17 



rt\pl,pl,pl) 



€=0 



€=0 



ri''\plplpl) 

A^''\pI,pI,pI) 



-8(n - 2) [(P1P2) + m^J (/'i + (n - 4)(pi - ps)^ + /«i,2) 
+2p3(pi - P2) Vi + 4(n - 2) (P1P2) + + Kx^2 - 2/«o,3) 

m2(pip2) 



+ (n-2) 

-n?{n - 2) 
0, 



P1P2 



+ 1 (P? - P2) ('^i.i - '^1,2) 

(P1P2)(P? +P2) 



- 2 



+ /«i,2 - 2k) 



g'r]CA 3 

^^t^Ca p\-pI [2(n- 1) 2 

P3 



€=0 (47r)"/2 32/C 

2(71-1 



/C 



(P1P2) 9^1+4 (P1P2) + 



(3.26) 
(3.27) 
(3.28) 



+ 



(P1P2) + (PlP3)'«l,l + (P2P3)/«1,2 + P3'«0,3 



+2p^(/7i + (n - 4) (ki,i + «;i_2) - 2 (ki^i + «;i_2 - 2^0,3) 
(n - 2) — (ki,i + Ki,2 - 2k) + (n - 4) (pi - P2) 



+{n - 2)m^ 
0, 



P1P2 

{PiP2){pI+pI) 



2 2 
Pi -P2 



P1P2 



ft^l.l - ^1,2 

Pi - pi 



(3.29) 
(3.30) 



q'^vCa 3 



{ps [(P1P2) +'m'^ 'Pi + 2}C(pi 



«=o (47r)"/2 4/C 

+ (PlP3)fi;i,l + (P2P3)«1,2 +P3'«0,3}- 



(3.31) 



In the limit of massless quarks, m — 0, the above results simplify considerably. First 



of all, in this limit, 



(1) 



s q 



1) _ ^(1) 



T^''' = = (this holds also in an arbitrary covariant 



gauge). Furthermore, cpo, cp2 (fo, Ki^i /to.i, ^2,1 ^0,1 and k ^ 0. 

3.3 Comparison with other papers 

As we have already mentioned, the contributions of the first diagram should coincide, 
up to an overall factor, with the one-loop contribution to the fermion-photon vertex in 
QED. Formally, to get the QED case from our expressions, we can put Ca =^ 0, Cp ^ 1, 
g e (the absolute value of the charge of the electron) . Therefore, our expressions also 
provide a one-loop correction to the QED vertex, in an arbitrary covariant gauge and 
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dimension. 

To get the expressions for the dimensionally-regulated four-dimensional case, we put 
n = 4 — 2e and expand our results (including the scalar integrals, see Appendix C) in 
e, keeping singular (l/e) and finite terms. In the Feynman gauge (^ = 0) and in four 
dimensions, we reproduce the well-known results from There are a few misprints 
in |[Tl| which were pointed out in |[12| (p. 1252). We agree with most of these corrections 
given m m, except for the following: Eqs. (3.12), (3.14) and (A19) of ^ are correct, 



they should not be changed. 

For an arbitrary value of ^ and in four dimensions, one-loop contributions to the QED 
vertex have been calculated by Kizilersii, Reenders and Pennington [O. We basically 



reproduce their results given by Eqs. (60)-(64), (67)-(74), and (87), except for and 
Te, Eqs. (69) and (72). Specifically, to get agreement with our results, we had to change 
the following^: in the sixth line of t^, 



m 



m 



)L'] 



and in the eleventh line of rg, 
3m2 



8A^ 



3m 
'8A2 



g'(/-F)[/(m'+A;')L+A;'(m'V)^'] 



For massless three-dimensional QED, results are given by Bashir, Kizilersii and Pen- 
nington in [|1^]. We reproduce their results for the longitudinal vertex, Eqs. (50) and 
(51), as well as the transverse parts,^ Eqs. (55), (57)-(60). 

Comparison with some other papers is given in Section 4, where the corresponding 
special limits are considered. 



3.4 Renormalization 



In the limit n 4 (e 0) the only function in the quark-gluon vertex which has an 
ultraviolet (UV) singularity at one loop is the A^^'* function. In an arbitrary covariant 
gauge, the UV- singular part of A^^'' reads 



(1,UV) 



(4vr) 



[1-0 {Cp-ICa) + 1{2-0Ca 



1 



+ 



{Airy 



1 

- + 

e 



(3.32) 



^There is also an obvious misprint in their (A14): the denominator \{k — 2)^ — m?] should read 
[[k — — m^]. We are grateful to the authors of for confirming all mentioned misprints. 

^However, we note that some factors of tt are inconsistent. The result for Jq (Appendix, Eq. (5)) 
should read (no tt). The right-hand sides of Eqs. (40) and (55) should have an extra tt. 

Finally, and K^^^ are related like Jo and J^*^-* . We are grateful to the authors of for confirming 
these misprints. 



19 



In the first line, the contributions from the first and the second diagram are explicitly 
separated. This result coincides with Eqs. (A.55)-(A.57) of Ref . ^O|0- 
The divergent parts of the two-point functions are as follows: 

.(^^-) = ^CHl-Og + ...), (3.33) 
/?(^-^^) = -j^m - f- + • • • V (3.34) 



(47r)2-^ 

The results for a and /3 agree with Eq. (2.5.138) of (our a and /? correspond to his 
(— -B) and mA, respectively), and with Eq. (11.55) of ||3^. Among the Xi functions, only 
Xo has a UV-singularity at one loop. 

Using our results we have checked that all other functions A^-^'' {i ^ 1), x^P (^ = 
1,2,3) and r/^'' have no UV-singularities. This was one of the important checks on self 
consistency of the calculation. 

To renormalize the above expressions, we need to subtract the 1/e poles, (possibly) 
getting some constant R instead, depending on the renormalization scheme: 

- + ...) ^ (i? +...). (3.37) 



In the MS scheme i? = 0, because (see Eq. (|2.7])) 

"l + 0{e^)\ , (3.38) 



7] = e '^'^ 



so that e '^^ and (47r)^ are absorbed by the MS re-definition of the coupling constant g^. 
This procedure can be also re-formulated in the language of the renormalization 



Z-factors, by analogy with Section VIII of [39 



4 Some special cases 

A few limits are of special interest: 

• the symmetric case, when pf = P2 = P3', 

^"^The contribution of the first diagram is also in agreement with Eq. (11.65) of p8| , whereas for the 
contribution of the second diagram there seems to be a misprint in Eq. (11.70) of | |38( |. Namely, in their 
notation, the factor (1 — ^) in Eq. (11.70) should read (1 + £,). Their ^ and e correspond to our (1 — ^) 
and 2s, respectively. 
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the on-shell limit Pi = P2 = (with or without the assumption that the vertex 
function is being sandwiched between Dirac spinors); 



• the zero-momentum hmit, when the gluon momentum vanishes (pa = 0). 

Since in all these cases we can put pf = P2 = p'^, we start by considering this as the first 
step towards all these limits. 

In the case pf = P2 = p^ some of the tensor structures ( |2.33| ) and (|2.36|) in the quark- 
gluon vertex become linearly dependent, namely: L2^^ and T2,^, Ls^^ and Ti^^, T^^^ and 
Tt^^. Moreover, according to Eq. ( |2.41| ), A4 and rg vanish. Therefore, the quark-gluon 
vertex in this limit can be written as [cf. Eqs. ( p.44| ) and (|2.45 )] 



Li^^Xi + L2,^A2 + L^,^i\z + T^,^ln + ^5,^X5 + Tt^^tj + Tg^^rg. (4.1) 



In fact, only the Li^^ contribution remains "non-transverse" in this limit. 



4.1 Symmetric case 

In Ref. 0] (see also in [^) the "symmetric" limit of the quark-gluon vertex, p\ = p\ = 
p'^=p'^ = —fi"^, has been examined (for the case of massless quarks, m = 0). The decom- 
position of the quark-gluon vertex is given in Eq. (2.34) of [Q. It basically corresponds to 
the naive decomposition presented in Eq. ( [A.l] ) of this paper. In Ref. ^ explicit results 
are given only for two scalar functions (out of twelve), 

Fi ^ hi and ri2 —h^ . (4-2) 

Taking into account Eq. (|2.41|) , in terms of A's and r's we get (see also in Appendix A): 



Ti = Ai + pgTs + {piP2)n 

Putting p\ = P2 = p\ = P^ (implying {P1P2) 
following one-loop results: 
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(4.3) 



-hp^) and m = 0, we arrive at the 



(1) 



g'rj (Cf - \Ca 



(47r)"/2 



-L 12 



(47r)"/2 

g^^ {Cp - \Ca) , 

(47r 
' (47r)"/2 



l"/2 



I [4(n - 4) - <] pVo + (3 - n + 
12 - 2(2n - 5)e + (n - 4)^'] /€o,3 - |(8 - <)pVo 
3 [{n -6) -in- 4)e] 



(4.4) 



^0 



Taking into account that the constant -R(l) used in can be identified as 
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n=4 



TC 



(4.5) 



(4.6) 



and expanding in e (keeping the divergent and finite in e terms), we arrive at the same 



nil 



In Eq. 



a2(^) = Im 



Li2(e 



is the Clausen 



result as given in Eqs. (2.36) of 
function. 

For the case of massive quarks (m 7^ 0), a similar "symmetric" limit has been con- 
sidered in Ref. 0, where the renormalization factor Zip was calculated at the one-loop 
order in the MOM scheme. For their calculation, only one of the scalar functions (namely, 
the one accompanying the 7^ matrix) was needed, Ti ^ hi, which is related to A's and 



r's via Eq. (|]3D. Putting 



-fi in our expressions, we arrive at the 



same result^as given in Eqs. (13)-(14) of 0. In particular, we have taken into account 
that the functions H and M introduced in Eq. (15) of are related to our functions ipi 
throughQ 



H 



2 / 2 2 2 
-P ^i{P ,P ,P 



n=A 



M 



2 / 2 2 2 
-P Mp ^P iP 



n=4 



(4.7) 



More details on these functions can be found in Appendix C. 



4.2 On-shell quarks 

The tensor structure of (f4.1|) does not change when we put = w?. Let us introduce 
the notation 

'^r,2 = Vi,2{m^^\pl) (4-8) 

for the functions in the on-shell limit. According to Eq. ( |C31D , in this limit the 
function (corresponding to diagram a) reduces to a two-point function and a tadpole. 



^2 



1 



\2{n - 3)fi:2,3 - {n - 2)k] 



(4.9) 



{n — A){Am? — p^) 

This is the reason why does not appear in the on-shell results for diagram a. 

The results for the relevant longitudinal and transverse functions in the limit p\ 
p\ = m? are presented below. 

4.2.1 Diagram a 



\(la)/ 2 



2 2\ 



g^7](CF-\CA) (n-2)(l-0 



(47r)"/2 



2(n-3) 



K 



(4.10) 



^^There is a misprint in Eq. (2.35) of Q: should read p^T 12 = —fi'^Ti2. Note that the e used 
in 0] has different sign, as compared to ours. 

^ Their a corresponds to our (f — ^), whereas their A denotes the ratio m,^/ii^ ~ —rn^/p^. Note that 
there is a misprint in the journal version of the resuh for the last contribution in Eq. (14), Cb'- in the 
very last term, In should read In j^j- In the preprint version of j^] this equation is correct. 



i+A- 



^■^ We also note a misprint in the large-mass ( 1/ A) expansion of the M function presented in Appendix C 
of Ref. §: the last available term, -2672/(11025A5), should read -1523/(6300A5). We are grateful to 
O.V. Tarasov for confirming this misprint. 
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(47r)"/2 (n-3)(n-4)(4m2 



X (n-3) {2-i){Am'-pi) + {n-h)ip, 



1^2,3 



-{n - 2)[Am^ - pI + 2{n - 5)^m^]K^ , 



(47r)"/2 (n-3)(n-4)(4m2-p2) 

(Cf — IC'a) 1 



(4.11) 

[(n-2)^-2(n-3)/^2,3] , (4.12) 



(47r)"/2 2(n-3)(n-4)pi(4m2 -p2) 
x| [(n - 2)(n - 3) - (n - 4)^]^^ [2(n - 3)«2,3 - (n - 2)k 
+{n - 3)(n - 4)(4m2 [2/,2,3 - (n - 2)k] | , 
(C^ _ ICa) 1 



(4.13) 



T^^''^(m^m^P3) = 



(47r)«/2 4(^ _ 3)m(4m2 - p^) 
+4(n - 3 - Orn^ [2(n - 3)^2,3 - (ri - 2)k] | , 
g'v {Cf - e 



n-2)e(4m'-p^)« 



(4.14) 



X 



T^^"^(m^m^P3) = - 
4.2.2 Diagram b 



(47r)'*/2 2(n-3)(n-4)m(4m2-p2)2 
|4(n - 5)m^ [2(n - 3)^2,3 - (n - 2)k\ - {n - 2){n - 4)(4to^ - pg)^! , 

(4.15) 

g^rj (Cf - ICa) 6 - n + (n - 4)^ 



(47r)"/2 (n-3)(n-4)(4m2-p2) 



[2(n-3)/^2,3-(n-2)^] (4.16) 



A^^''^(m^m^P3) 

A^^^^(m^m^P3) 



^'r^C^ 2-e 



(47r)"/2 4(n-3) 

1 



(47r)"/2 16(n-3)(4m2-pi)2 



[{n - 3)ko,3 + (ri - 2)k] , 

8(n-l)(n-3)(2-Om^pl(^r 



(4.17) 



-2(n-3)^e P3[2(n-3)m^+p,^]^f - 2(n-3) V[4(n-4)m^ + 3p2]fi:o 3 
+4(n - 3)(2 - 0[4(n - 2)m2 +p2]^o,3 + (ri - l)(n - 2)(n - 3)eV3^ 
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(U) / 2 2 2\ 



/ 2 2 2\ 



(Ife)/ 2 2 2\ 



(4.18) 



g'yCA 1 

(47r)"/2 32(n - 3)m(4m2 - p|) 



|2(n-3)(n-4)fp^^4m2_^2)^os 



+ 



(4.19) 



4(n - 1) + 2(n - 4)^ - {n - 2){n - 3)^^ 

X 4(n - 3)mV3V5i' + 8(?2 - 3)m^/to,3 - (n - 2)^3/? 

+4(n - 3)^^^(4m2 - ^3)^^0,3 - 4(n - l)(n - 2)(4m2 -p^)/^ 

^^^^A 1 

(4vr)n/2 I6(n - 3)j9i(4m2 - p2) 

X |2(n - 3)^(4 - 0(4m2 - pj) - (n - 4)/^o,3 

- 2(2-0 - (^-3)^2] [4(n-3)mVr + 2(n-3)Ko,3 - (n-2)/5 (4.20) 

q'vCa 1 

(47r)"/2 8(n - 3)m(4m2 - p^) 

+^[2 + {n- 3)^]m^ [4(n - 3)mVr + 2(n - 3)ko,3 - (n - 2)k 
-2{n - 3)[6 + 2(n - 5)^ + ^ V(4m2 - p2)^os| ^ ^^ 21) 

^^^^A e 

(4vr)"/2 I6(n - 3)m(4m2 - ply 
X |2(n - 3)(4m2 - pi) [im^ + 2{n - 5)^m^ - {n - 4)^pl\ y?"' 
+ (n-l)[2 + {n- 3)^] 4{n-3)mY3VT + 8{n-3)m^Ko^3 - {n — 2)p^K,\ 
-4(n - 3f^{Am^ - ^3)^0,3 - (n - 2) [6 + (n - 3)^ (4m^ - pl)R^ , (4.22) 
1 



■|(n - 2)(n - 4)04m^ - P3)k 



(4.)"/. 8(„ - 3)(4,n= - Pi) l^<" - - «><^'"^ - 

12 + 2(2n-7)^ - (ra-3)^^] [4(n-3)mVr + 2(n-3)Ko,3 - {n-2)R 

(4.23) 



When P3 — s> 0, the only problem which may arise in Eqs. ( [4.1(]| )- (|4.23| ) concerns the 
functions rj^"'' and t^^''^ containing p| in the denominator. However, the coefficients 
of {pD^ in the numerators of Tg^""* and rj^''^ are proportional to [2/t2,3 — (n — 2)k\ and 
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Ko,3, respectively. If we take into account that K2,3|p2^o ~ li''^ ~ '^)'^ ^^0,31^2^0 — 
(massless tadpole) , we see that the hmit p| — > is regular for r^^^ . 

4.2.3 Dirac and Pauli form factors 

Let us consider again the on-shell limit pi — P2 — without putting any further con- 
ditions on P3. If we recall that the "physical" quark-gluon vertex should be sandwiched 
between physical states obeying the Dirac equation, 

It(-pi)r^it(p2), with It(-pi)^i = -mir(-pi), ^2u{p2) ^mu{p2), (4.24) 

then (using the above Dirac conditions) we arrive at the standard decomposition 

u{-Pi)'^Mp2)=Fi{pI) u{-pi)j^u{p2) - ^^2(^3) u{-pi)a^yplu{p2) 



U{-Pl)lt^u{p2) + -^F2{pl)u{-Pi){pi-P2)^u{p2), 



(4.25) 



where Fiij)^ and -^2(^3) are often associated with the Dirac and Pauli form factors, 
respectively. In terms of the (modified) A and r functions we get 



Fi + F2 = Ai +p^T3 - 2mT5 + |(p1 - 4m2)T8, 



1 



F2 = -2m\2 + A3 - T5 + IP3T7 - mrg. 



2m 

Using our results we get 

g'^r^ (Cf - \Ca 



(4.26) 
(4.27) 



(Fi + F2)(i") = 



«2,3 



(4.28) 



(47r) {n - 3) (n - 4) {Am? - pi) 
X {(n - 3) [2(n - 3)p^ - (n^ - 9n + 22) (4m^ - pi) 

-2{n - 2) [2(n - 5)m^ + pi] k} , 

^Uin - 3)p^(2mVr + ^os) 

(47r)"/2 8(n-3)(4m2-p2)\ V '^^^ ^1 

+{n - 3) [S + 2(4n - 13)^ - {n - 4)^] {^m^ - P^)«o,3 

-8(n - 2) (5m2 - p|)k| , (4.29) 



2m ^ 



(Cf - ICa) (n - 5)m 



(47r)"/2 



(n — 3)(4m2 — pf) 



[2{n - 3)k2,3 - (n - 2)^] , (4.30) 
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X 



.n 



3)pI 12m^ + (n - A)pI y^f + 2(n - 3) 



(n 



+ - 2) {n- 6)(W - pi) - {n - l)p\ 



K 



(4.31) 



The results for the (QED-hke) diagram a agree with Eqs. (5.69) and (5.70) of ^ . 

In the massless hmit, m — 0, and in Feynman gauge, ^ = 0, our results can be 
compared with those of where renormalized results for the quark self energy, S, and 
the quark-gluon vertex functions, and (their notation), are collected in Tables B.I 
and B.II, respectively. The vertex functions are given for off-shell gluon {pf = P2 = 0, 
p| 7^ 0) and the out-going quark being off-shell [pf 0, p^ = p^ = 0). We confirm these 
results. 



4.3 Zero- momentum limit 

Let us consider the zero-momentum limit p3 = 0, with p2 = —pi = p being off shell. 

We can proceed in two steps. The first step, putting pf = P2 = (without putting 
p^ = m^), has been already done in Eq. ( [4.11) . As the second step, we should put pi = 
(which implies ps = 0). In this limit, the quark-gluon vertex becomes 

r^|p2=_p^=p^ P3=o = ^i.M^i + ^2,mA2 + Ls^i^h = 7mAi + 4:p^^X2 - 2pf,X3. (4.32) 

One can check that the corresponding scalar functions are regular. 
The substitutions for the relevant integrals are^"^ 

^2(1, 1, l)U=-,,.p,,3=o = -^^2 [(^ - 3)>/2(0, 1, 1) - Ml, 1, 0)] , (4.33) 

n — 2 

MlAM,,=-,,.,,,,=o = ^ ^2(0,1,0) , (4.34) 

^1(1, 1, l)U=-...p,P3=o = J^^^2 {(^ - 2)^i(0, 0, 1) 

-(n-3)(p2 + m2)Ji(0,l,l)} , (4.35) 

MlAM,„^,,,,=o = 0- (4-36) 
The scalar functions from ( [4.32| ) in this limit are 



^"^In general, one should also make sure that the next term of the expansion in p§ does not contribute, 
which may happen when appears in the denominators. In this calculation, we did not need such p| 
terms. 
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g^T] (Cf - ^Ca) (n-2)(l -0 



+ 



(4.38) 



[Cf - ^Ca) {n - 



(47r)"/2 2p2(p2 _ ^2) 

X I (n - 4)p^ + (n - 2)m^ k,i{p^) - {n - 2)m^Kj 

q^vCa 2 - e 



(4.39) 



n — 2)(m^ — 3p^)m^K 



+ 



4(n - 3)mV + (^ - 2)(p^ - m^)(3p^ + m^)] fti(p^)} , (4.40) 



Af^(p^P^O) 



(^^ryC^ (n-2)(2-0 



{ [(n - 4)p2 + 



nm 



(47r)"/2 32(p2)2(p2 _^2) 

- 4(r2 - 3)m^ + (n - A){p'^ - m^){p'^ + 3m^)] Ki(p^)} , (4.41) 

2 2 n^ 9'^VCa m{n-l) r r 2 , / 2l ( 2\ 

^3 >0) = (^^^ 4p2(p2_^2) [[{n-2)m + (n - 4)p J /ti(p ) 

-(n-2)m2K} . (4.42) 

If we now consider the on-shell case, i.e. put = m? , we need to be careful, since 
the above expressions contain {p^ — m?) in their denominators. Using two terms of the 
expansion of Ki(p^) in 5 = {m? — p'^)/m?', 



n-2 
2(n-3)' 



-Kim 



1 + U + 0{6^) 



(4.43) 



we see that the pole at = (5 = 0) is canceled. In this way, we arrive at 



(47r)"/2 2(n-3) 
g'^f] in — 2)fi; 



C^-lC^) (l-0 + iC^(2-0 



(47r)"/2 4(n-3)m2 
(yf2r^ (tt, — 2)/5 



(47r)"/2 4(n-3)m 



(4.44) 

Cf - \Ca) (1 - + iC'^lS - 0] , (4.45) 
Cf - \Ca) (ri - + IC'a (ri - 1)1 • (4.46) 



These results coincide with those obtained by first taking the on-shell limit p\ = p^ 



m 



(see Eqs. ( CTD -( gl2D and ( gTzD -f gloD above) and then putting p^ = 0. If we 
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assume that the vertex ( ^.32| ) is sandwiched between physical states ( [4.24| ), we see that 
both structures containing can be transformed into 7^, so that the "effective" vertex 
becomes 

7^[Ai(m^ 0) + 4m^A2(m^ 0) - 2mA3(m^ 0)]. (4.47) 
Substituting our results ( [4.44|) -( p:.46D we obtain for the one-loop contribution to Eq. ( [4.47|) 



7m 



g'^T] (n- l)(n-2) ^ _^ 



(47r)"/2 2{n - 3) 



(An] 



2-s 



^(e)|3||cHm^)-^ 



(4.48) 



where we have taken into account the definitions of rj ( ^.7[ ) and k ( ^.151) . We see that 
the result ( [4.48|) is gauge independent, and it does not contain Ca- For the QED case, 
it coincides with Eq. (22) of . 

In the massless (m = 0) case, we can compare our results for the zero-momentum 
limit with those given in The definition of the zero-gluon-momentum vertex is given 
in their Eq. (A4) (the lower equation). It is proportional to (their q our p) 



q2 



(4.49) 



If we consider our results in the off-shell zero-momentum limit, we see that on putting 
m = the function A3 vanishes. The remaining two functions can be mapped into r3 4 

as 

Af'") ^ r3 + r4, -4p2A2 ^ r4, (4.50) 



or, vice versa, 



-4p'A2, 



(4.51) 



where the superscript "(ren)" means "renormalized" . 

In the massless case, taking into account that the massless tadpole function k 
vanishes, we get 



m=0 



(1) 



g^ 7] n — 2 



(47r)"/2 



(4.52) 



4p^A 



2T(i) _ 9^V (n-2)(n-4) 



(47r)"/2 



{P') [{Cf - ICa) (1 - + \Ca{2 - 0] .(4.53) 



Putting p2 = /i^ (as in ||^) means that the logarithms ln{p'^ / fi'^) (in our case, ln{p'^), 
since we imply /i = 1) vanish. Therefore, we should substitute kq =^ 1/e + 2 (see, for 
example, Eq. (2.15) of P^)- Omitting the 0{e) terms we see that 



(1) 



4/A?) 



167r2 



Cp — \Ca 



1-0 + !Ca(2-0 



167r2 



2(Cf-\Ca){1-0 + \Ca{2-0 



(4.54) 
(4.55) 
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where we have exphcitly separated the contributions of diagrams a and b. To perform 
the renormahzation (in the MS-scheme), we need just to subtract the 1/e term, i.e. 



(l,ren) 



167r2 



- he, 







(4.56) 



Now we can compare our Eqs. ( [4.55| )-( p:.56| ) with one-loop results for r3^4 presented 
in Eq. (A12) of p. Taking into account that their ^ corresponds to our (1 — ^), and also 
that their a = / (47r), we see that we are in agreement with Moreover, their result 
for the quark self-energy function A{q^) (see Eq. (A9) of 0) is in agreement with our 
a{p^). Note that /5(p^) = in the massless case. 



5 Conclusions 

In this article we have given results for the one-loop quark-gluon vertex in an arbitrary 
covariant gauge and in an arbitrary space-time dimension. The calculation was carried 
out with massive quarks. 

To calculate the quark-gluon vertex, we have decomposed it into longitudinal ( p.32| ) 
and transverse ( f^.35D parts, V^^^ and V'^^'^ (like the decomposition in QED |TT[). Alto- 



gether twelve scalar functions (four A's and eight r's) are needed to define the quark-gluon 
vertex. We found that the function A4, the coefficient of cr^jy(pi —P2Y^ which is absent in 
QED, does not vanish in QCD (see Eq. ( p.l2| )) and contributes to the non-Abelian part 



of the Ward-Slavnov-Taylor identity ( p.29|) . The general results for the longitudinal 
parts (A's) are given in Eqs. (|3l5|)- (|3.12|) (arbitrary gauge), and results for the transverse 



parts (r's) are given in Eqs. (|3.16|) - (|3.31|) (Feynman gauge) and in Appendix E (arbi- 



trary covariant gauge). Using recurrence relations (see Appendix B), all results have 
been expressed in terms of integrals with powers of propagators equal to zero or one. 
Only two non-trivial scalar functions are required, ^iiPiiP^iPD-, where i = 1 or 2 counts 
the number of massive propagators involved (see Section 2.1 and Appendix C). 

Starting from the general off-shell expressions {p\ ^ m^, p\ 7^ m^, p\ 7^ 0) in an 
arbitrary space-time dimension, ra, for the longitudinal and transverse parts of the vertex, 
we have derived results for the on-shell limit {p\ = p\ = m^) which are also valid for 
an arbitrary n (Section 4). Note that if we started from the off-shell results expanded 



around n = A (similar to the results of 12| for the QED case), we would get infrared 



divergences from the on-shell-divergent logarithms. Keeping the arbitrary space-time 
dimension, we see that the corresponding infrared divergences result in extra poles in 
e = (4-n)/2. 

Various special cases of the general results were compared with those of Refs. 
|,|n|,|g,|T§|g (for details, see Sections 3 and 4). 

At the one-loop level, the Ward-Slavnov-Taylor identity for the quark-gluon vertex 
can be split in an Abelian and a non-Abelian part, Eqs. ( |2.28D and ( |2.29| ). The Abelian 
part is similar to the Ward-Fradkin-Takahashi identity in QED |3^, whereas the non- 
Abelian part has a nontrivial contribution involving the quark-quark-ghost-ghost vertex 
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( p.25| ), which can be described by scalar functions Xi = 0, ...,3). One-loop results for 



these functions are presented in Appendix D. Using the results for these x functions, 
and those for quark and ghost self energies, we have checked that our results for the 
longitudinal parts of the vertex (Ai, A2, A3, and A4) satisfy the WST identity for arbitrary 
n and ^, as they should. 

In principle, some techniques which can be used for the calculation of the two-loop 
off-shell quark-gluon vertex, at least in the m = case, are already available [p^ , 
although the problem of higher powers of irreducible numerators is still difficult for 
algorithmization. For special limits, the calculation is very similar to the three-gluon 
vertex, which was calculated at two loops in (the zero- momentum limit) and in 
(the on-shell case). 

Note. The results of this paper are available in REDUCE format at 
tittp : //wwwthep . physik . uni-mainz . de/Publications/progdata/mzth9963/| 
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Appendix A: Decomposition of the vertex 

The general quark-gluon vertex can be expressed as (see |TT| , p!^ ,p5|) 



+pi,,hh% + piMi^i/ig + Piiiixh^Q + piiiiih^w + iii'i>\h^\2-, (A.i) 

where hi = hi (pj ,pI,pI). 

The longitudinal and transverse functions of Eqs. ( |2.32| ) and ( [^.35| ) are related to this 
naive basis as follows: 

Xi = hi - \{p2Pz){hQ + hj) - \{piPz){h^ + /ig) + \{pI - 2pl)hi2, 
A2 = 2(p2 \ p2-^ [iP2P3){h6 - h-r) + {piP3){hs - hg) - plhu] , 

A3 = 2~^{(P2P3) [h2 + hi + {piP2)hio] + (piPs) [/l3 + /^s + (P1P2) /^ll] } , 

Pi ~ P2 



A4 = I [-/i4 + /i5 + {P2P3)hio + {piP3)hu] , (A. 2) 
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Tl = 2 ^ 2 [^2 + /is + ^4 + /^S + {PlP2){hiQ + /in)] , 
Pi ~ P2 

'^2 = TZT^ 2t(^6 -h7 + h8-hg + 2/ii2), 

^3 = -i(^6 + /l7 + + /ig), 
r4 = 2 ^ 2 (^10 + ^11), 

Pi -P2 

T6 = i(^6 + hY-hs-hg- 2/li2), 
= 2 ^ 2 [(P2P3)/ilO + (PlP3)/ill] , 

Pi -pi 

Ts = -/il2, (A.3) 

where = Xi{pl,pl,pl), = n{pi,P%pl), and hi = hi{pl,pl,pl). 
Inverting these relations, one finds: 

hi = Xi+ TspI + Tq{pI - pI) + Ts{piP2), 

h2 = -A3 + A4 - n{pips) + r4(piP2)(PlP3) + + T7{pI - 
/l3 = A3 - A4 + [ri - n{piP2)] {P2P3) + Ts - T-jipl - 1^3), 

hi = -A4 - Ts + ir7(p^ - pI), 

/l5 = A4 - Tg - ir7(p? -p^), 

he = X2- T2{piP'i) - r3 + re, 

/i-7 = -A2 + T2{piP^) - T-3 + Te - Ts, 
/l8 = -A2 + T2{P2P^) - ^3 - Te + Ts, 

/ig = A2 - r2{p2Pz) -Ts-Tq, 
hio = -nipips) + tt, 

hu = nip2P3) - -Ty, 

hi2 = -Ts- (A.4) 

Appendix B: Recurrence relations for scalar integrals 

To calculate scalar integrals with higher (integer) powers of denominators, a recurrence 
procedure based on the integration-by-parts technique ||2^ has been used. For the three- 
point integrals, one can follow an approach described in (see, in particular. Section 4 
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of where the massive case is discussed) . Using the recurrence relations, we can reduce 
all integrals to the master integrals Jj(l, 1, 1) {i = 1,2) and a few two-point integrals. 
All basic integrals are discussed in Appendix C. 



B.l Two-point integrals 

All two-point integrals occurring in this paper can be associated with certain special 
cases of three-point integrals J2 and Ji (defined in Eqs. ( |2.4| ) and (p.5|), respectively), 
when one of the indices z/j vanishes. The masses of internal particles (mi,m2) can be 
equal to (m, m), (m, 0), (0,m), or (0,0). 

For two-point integrals with arbitrary masses, 

J(^)(.„.2|m„m2) = /d"g[(^_^)2_^2p(^2_^2).. ' (B-1) 



such procedure has been described in detail in Appendix A of For positive Ui, it is 
enough to apply 

J^^^(z/i + 1, z/2|mi,m2) = --^|[(n-z/i-2z/2)(p^-m^) + {n-3ui)ml]j'''^\ui, i'2\mi,m2) 

-uiip^ - m\- ml)J^'^\ui + 1,U2 - l|mi,m2) 
-2i^2mlJ^^\ui - l,z/2 + l|mi,m2)j, (B.2) 



together with a similar equation with (z/i,mi) ^ {i^2j^2)- In Eq- 

A = A(m^, ml, p^) = 4m^m^ - (p^ - mj - mlf . (B.3) 

In our calculation, A may be equal to p'^{4:m? — p"^), —{p^— m?Y "{p^Y- 

Note that the sum of indices Vi for any of the integrals on the r.h.s. of Eq. (|B.2| ) is less 
(by one) than such a sum for the integral on the l.h.s. Therefore, J'-^-* with any (positive) 
integers ui and z/2 can be reduced to J^^^ with vi = 1/2 = 1 and tadpole integrals. 
The reduction of tadpole integrals is trivial: 

Ji(0,0,Ar)= ^^^^^^ (1-1)^,1 (-"^')'-^ ^1(0,0,1) (B.4) 

(and similarly for J2(A^, 0, 0) = J2(0, N, 0)), where (a)-,- = r(a + j)/r(a) is the Pochham- 
mer symbol. 

B.2 Three-point integrals J2 

For the integrals J2 with positive Vi (with at least one of them being greater than one), 
the following solution of recurrence relations can be used: 



J2(.i,.2,.3 + 1) = ^ 



{yi-yz)plpl + i.^2-yz)plpl + 2(n-i/i-i/2-2z/3)p3(pip2) 
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+2(z/2 - z/i)m - P2) - {2n - Ui - U2 - Qu^)m 



J2(z^i, 1^2, 1^3) 



-pUpI - 4m^)[z^i ^2(^^1 + 1, 1^2, ^^3-1) + l^2M'^li ^^2 + 1, 1^3-1)] 



X [Z/2^2(Z^1 - 1, Z^2 + 1, 1^3) + l'3J2{l'l " 1, ^^2, ^^3 + 1)] 



(B.5) 



J2{Vl + 1,1^2,1^3) 



1 



2(n - z/i - z/2 - 2z/3) plip2P3) - m^iPiPs) ■M'^i, '^2, 1^3) 



2U1M2 

+ {pI - m^) (z/3 - Z/i)p^ + (z/2 - Z/3)P2 + (i^l - 1^2)"^^ -^2(Z^1, ^^2, T^s) 



+ 



P3i.P1 - m)-2m {p^ -P2) 



x[z^lJ2(«^l + 1, Z^2, Z/3 - 1) + Z/2J2(Z^1, Z^2 + 1, Z/3 - 1)] 

+ (Pl-m^)(p2-"^^)[^^l^2(^^l + l, «^2-l, ^^3) + 1^3^^!, ^2^1, «^3 + l)] 

-(p?-m2)2[zy2J2(z/i-l,Z/2 + l,Z/3) + Z/3J2(z/i-l,Z/2,Z/3 + l)]j , (B.6) 



and also an equation for J2(^i, ^^2 + 1, ^'3) which can be obtained from ( p.6|) via (z/i,p^) 
(^2,^2)- The quantity AI2 is defined in Eq. (|3.3| ). 

In fact , for our calculation we needed only the ui = 1/2 = 1^3 = 1 case of Eq. ([B.5|) , 



^2(1,1,2) 



2(71 -A)pi ipiP2)-m' J2(l,l,l)-2p^(p^- W)J2(2,1,0) 



2A^2 

(pI - mV3 + '^{pI - pIW] [-^2(2, 0, 1) + J2(l, 0, 2)] 



{pi - m')pi - 2{pi - pi)m' [J2(0, 2, 1) + J2(0, 1, 2)] 



(B.7) 



where we have taken into account the obvious symmetry J2(z^i, ^2, 0) = ^2(^^25 J^i, 0). We 
note a remarkable fact that the coefficient of 72(1; 1) is proportional to (n — 4). Since, 
in the off-shell case, J2(l, 1, 1) is finite in four dimensions, this means that we do not 
get any nontrivial function in finite parts of all triangle integrals with higher z/j. In fact, 
this property is valid for arbitrary internal masses |^ (a special case has been discussed 
in il). 
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B.3 Three-point integrals Ji 

For the integrals Ji with positive Vi (with at least one of them being greater than one), 
the following solution of recurrence relations can be used: 



Jl(l/1, 1/2, f-i + 1) 



Jl(z/i, Z/2, Z/3) 



Jl(z/i + 1, z/2, z/3) 



2z/3A^i 

+ (2r2 — 3z/i — 3z/2 — 2Z/3) 
+ (^2 - m^)[z/i Ji(z/i + 1, z/2 - 1, z/3) + z/3 Ji(z/i, z/2 - 1, z/3 + 1)] 

+ - m^)[z/2Jl(z/i - 1, z/2 + 1, z/3) + Z/3Ji(z/i - 1, z/2, Z/3 + 1)] 
-P3[z/i Ji(z/i + 1, z/2, z/3 - 1) + Z/2Ji(z/i, z/2 + 1, Z/3 - 1)]| , (B.8) 

2(n - Z/i - z/2 - z/3) Pi(p2P3) +"^^(PlP3) Jl{Vl,V2,V-i) 



-{VlPl - Z/2P2 - l^3P3)PlJli^l^ ^2, 1^3) 

2uipI + 2z/2(piP2) + (z^2 - ^^3)^3 " (^1 ^ ^^2)^^ Ji(z/i, z/2, z/3) 

+ +P3m^)[z/iJi(z/i + l,Z/2 - l,Z/3) + Z/3Ji(z/i,Z/2 - 1, Z/3 + 1)] 

+pI{pI - m^)[viJi{vi + 1, z/2, z/3 - 1) + z/2Ji(z/i, z/2 + 1, z/3 - 1)] 

-(p?-m2)2[zy2Ji(z/i-l,Z/2 + l,Z/3) +Z/3Jl(z/i-l,Z/2,Z/3+l)]| , (B.9) 



and also an equation for Ji(z/i, z/2 + 1, ^^3) which can be obtained from ( [B.9| ) via (z/i,p^) ^ 
(z/2,P2)- The quantity A^i is defined in Eq. (|3.2| ). 

The "highest" integral Ji which occurred in our calculation was Ji(2, 2, 1). Using 
Eq. ( p.9| ) and the (z/i,Pi) ^ (^^2,^2) symmetry, we get 



Ml + m'p^ [Ji(2, 1, 1) + Ji(l, 2, 1) + Ji(l, 1, 2)] 



Ji(2,2,l) 



2piM, 



+ {n-b) \fi{piP3) + m^{P2P3)\ Ji(2, 1, 1) + (n-5) [p\{p2Pz) + m^{piPci)\ Ji(l,2,l) 
-{pl-rrt'f [Ji(0, 3, 1) + | Ji(0, 2, 2)] - {pl-m^f [ji(3, 0, 1) + | Ji(2, 0, 2)^ 

+P3(p? + P2-2m2) [j^(3,i,o) + iJi(2,2,0)] } . (B.IO) 

For the integrals Ji(l, 1, 2), Ji(2, 1, 1) and Ji(l, 2, 1), direct application of Eqs. (p.8|) - 
;9D yields: 



Ji(l,l,2) 



2Mi 



2{n-A) {piP2)+m^ Ji(l, 1, 1) - 2^^ 7^(2, 1, 0) 
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+ (p^-m2)[Ji(2,0,l) + Ji(l,0,2)] + (p?-r/i2)[Ji(0,2,l) + Ji(0,l,2)]|,(B.ll) 

Ji(2, 1, 1) = ^;^{2(n - 4) [pI{p2P3) + m\p^p^)] Ji(l, 1, 1) + 2pl{pl - m^)J,{2, 1, 0) 

+{Mi+my,) [Ji(2, 0, 1) + Ji(l, 0, 2)] - {pj-m'f [MO, 2, 1) + Ji(0, 1, 2)] | 

(B.12) 

and similarly for Ji(l, 2, 1). 

B.4 Integrals with numerators 

We also need some integrals with negative powers of denominators (i.e., when the corre- 
sponding denominator is in the numerator). Such formulae can be obtained in a standard 
way, via representing the numerators in terms of the external invariants. For the integrals 
J2, we only needed the cases when one of the z/j was equal to (—1), 

M^i^ ^2, -1) = - (P1P2)] M'^i^ ^2, 0) + J2{i^i, 1^2 - 1, 0), 

J2(z/i, -1, Z/g) = -{pI)~^{P2P3) [M^U 0, Z/3 - 1) - ^^2(2/1, 0, Z/3)] - (P1P3) ^2 , 0, U^), 
-^2(-l, J^2, i^s) = -{ply^iPlPs) [-^2(0, U2, i/3 - 1) - J2(0, U2, i/3)] " (P2P3)-^2(0, U2, ^3). 

We also hst analogous results for Ji: 

i^2, -1) = -(P3)~^{P3[(P1P2) + m^]Ji{l/i, 1/2, 0) + (piP3)^l(z^l - 1, 1^2, 0) 

+ {P2P3)M^1^^2 - 1,0)| , 
Jl(z/i, -1,1/3) = -(P2)"^{b2(PlP3) +"^^(P2P3)]^l(i^l, 0,1/3) 

+ (P2P3)^l(i^l, 0, i/3 - 1) + (piP2)^l(j^l - 1, 0, i/3)|, 
Jl(-1, 1/2,1/3) = -(p?)~^|[Pl(P2P3) +m^(piP3)]^l(0,I/2,i/3) 

+ (PlP3)-^l(0, i/2, i/3 — 1) + (piP2)-^l(0, i/2 — 1, 

Moreover, an integral with i/3 = —2 has occurred, 

Ji(i/i,i/2,-2) = |[(piP2)+m']'-(n-l)-^/c|ji(i/i,i/2,0) 
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+2(p2)-i|(pip3)[(piP2) + (n- l)-i/c|ji(z/i - l,z/2,0) 

+2(p2)-i|(p2P3)[(piP2) + m^] + (n - l)-i/c| Ji(z/i, - 1, 0) 



2\-2 



+ ^3) 



(PlPa)' - - l)-'/C 
{V2Vzf -{n- l)-i/C 



Ji(z/i -2,z/2,0) 



Ji(z/i,z/2-2,0) 



with /C defined in Eq. ( p.ll) . 



(PlP3)(P2P3) + - 1) 



Jl(z/l-l,Z/2-l,0) 



Appendix C: Basic scalar integrals 

As discussed in Section 2.1 (see also Appendix B), all results for the scalar functions 
(Aj, Tj, etc.) can be expressed in terms of three-point integrals Ji(l, 1, 1) and 72(1, 1, 1), 
two-point integrals ( |2.1CI| )- (p.l3|) and a tadpole integral ( |2.14]) . Remember that we are 
interested in results which are valid for arbitrary values of the space-time dimension n. 



C.l Two-point integrals 

The massless two-point integral is trivial, 



Ji(l,l,0) = Jo(l,l,0) = ivr'^/2 {-Vlr'''-^ 



(f - 1) r (2 - f ) 



Yin - 21 



(C.l) 



Two-point integrals with one or two massive lines can be expressed in terms of the Gauss 
hypergeometric function (see, e.g., in ||2l 



Ji(0,l,l) = J2(0,l,l) = 2iWMm2)"/^-2^^fi-il 2F,f^' 

n — 2 \ n/Z 



J2(l, 1, 0) = ivr"/^ (m^)"/2-2 r (2 - f ) 2^1 



1, 2-n/2 
3/2 



Am? 



, (C.2) 



(C.3) 



The results for the integrals Ji(l, 0, 1) = 72(1, 0, 1) can be obtained from ( |U.2| ) by sub- 
stituting pI — > P2- 

In three dimensions, Eqs. ( |C.2| ) and (|C.3|) yield (cf., e.g., in [^) 

Jl(0,l,l)U3= J2(0,l,l)U3 = ^/f4l , 



1 ^1 \n=3 



m \ Am? 



(C.4) 
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with 



7^ arctan \/^^ , 2; < . 



(C.5) 



Around four dimensions (n = 4 — 2e), they are singular (see, e.g., 



Ji(0, 1, 1) = J2(0, 1, 1) = i7r'-=m-'^r(l + £)<!- + 2 + 



1 



J2(l, 1, 0) = m'-'m-'' V{l + e)l- + 2-2f 



11 11 
m Pi m —pi 



Pi 



In- 



0{s) 



pI 



0{e) , 



(C.6) 
(C.7) 



— ^m? 

with the same function / as in Eq. ( p.5|) . Above the corresponding threshold {p\ > m? 
or p\ > 4m^), these functions acquire imaginary part, whose sign is defined by the causal 
prescription + iO (for details, see Appendix A of ||46[|). 

In the limit pf = P2 = P^ (in particular, in the on-shell case), the combinations 

J2(l,0,l)- J2(0,l,l) 
P2-P1 



should be treated in the following way: 

rj2(i,o,i)-j2(o, 1,1)1 



pI-pI 



p\=pI=p'^ 



Using Eq. (|C2|), and also the fact that 



^{^2(0,1,1)1,..,.} . 



az 



a, h 



ah 

— 2-ri 
C 



a+ 1, 6+1 
c+ 1 



(C.9) 



we obtain 



j2(i,o,i)- j2(o,i,i: 

pI-pI 



In the on-shell case, we get 

^2(0,1,1) 



n-2 



2rn?{n — 3) 



^2(0,1,1)1,.=,.. (C.IO) 



^2(0,1,0). (C.ll) 



Moreover, for some functions we need the expansion of this integral in 5i = {m? —p\)/m? 
up to the linear term. 



^2(0,1,1) 



n-2 

2iv?{n — 3) 



J2(0,l,0) \l + \5i + 0{5l 



(C.12) 
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C.2 Three-point integrals 



For the triangle integrals Ji and J2, results in arbitrary dimension, and also for any 
powers of the propagators, can be presented in terms of multiple hypergeometric func- 
tions [^. For our purposes, we need integer powers of propagators. Using recurrence 
relations based on the integration-by-parts technique ||2^ , all scalar integrals can be 
reduced to Jj(l, 1, 1) and two-point functions. 



C.2.1 General off-shell case 



Transforming Feynman parametric integrals (for example, using the Cheng- Wu theorem 
we can present Ji(l, 1, 1) and J2(l, 1, 1) as 



Ji(l,l,l) = -i7r"/2r(3- 



d^ dr] 



J2(l,l,l) = -i7r"/2r(3-f) 



{ { (l+^-hr^)" ^[m'^r]{l+^+'r])-r]pl-^'r]pl-^p] 



d^ dr] 



2i3-n/2 



(C.13) 







(C.14) 

In the three-dimensional case {n = 3) the denominator {I + C.+tj) disappears in both 
integrals, and one can easily integrate over ^, and then over rj. This yields 



^i(i,i,i)L= 



ITT 



P3M1 



In 



m + JMi/pI 



J2(l,l,l 



ITT 



5 -^5 ^J\n=3 



In 



m 

m{2m? —p\— pI) + 
m{2m? — Pi — P2) — \J M.2 



m,/pI^ 



(C.15) 



(C.16) 



where A^i and J^2 are defined in Eqs. ( |3.2| ) and (|3.3| ), respectively. If M.i/p\ or M.2 
is negative (in Ji and J2, respectively), the logarithms should be substituted by arctan 
functions (see Eq. ( p.5| )), which correspond to the hmiting cases of the result obtained 
in (see also in ^3|) for three-dimensional three-point function with arbitrary masses. 
Note that in the massless case we arrive at the well-known result |^ (see also in ||14| ) 



-^1,2(1, 1,1 



1 "-1 ^J\n=3, m=0 



Jo(ia,i 



ivr 



, \n=3 



111 

-vivivi 



(C.17) 



In four dimensions (n = 4), we can also integrate over ^. Then, performing the 
remaining 7] integral, we arrive at the known results in terms of dilogarithms |^ . 
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C.2.2 Symmetric case 

The "symmetric" case is of a certain interest (see e.g. in 0), when all external invariants 
are equal, pf = pi = pI = — A*^- Then we obtain 



■^2(1, 1, l)l„=4, p2 = -^t2 



ITT 



drj 



fi^ J 1 + 1] + T]^ 



In 



ivr 



d?7 



/i^ J 1 + 7] + 7]-^ 



In 



(m2 + + r/)2 



m2(l + r])^ + /i^?] 



(C.18) 



(C.19) 



These parametric representations are equivalent to those for the H and M functions given 
in 0. Note that n^Ji itt^H and /i^J2 ivr^M. These integrals can be evaluated in 
terms of Clausen functions as 

,2 



ivr 



+012 (vr - 29,^) | 



+ 2CI2 - + 20^1 ]+Ch{-- 26, 



TT 



7sl 



(C.20) 



-^2(1, 1, 1)L=4, p2 = _^2 



2i7r^ 



CI2 



TT 



7r 



where 



tan^si 



tan ft 



s2 



3/i2 



(C.21) 



(C.22) 



This gives analytical results for the H and M functions from In the massless limit 
(m —>■ 0), 6si = Os2 = 7i"/6, and, remembering that CI2 = |Cl2 (|), we reproduce 



the well-known result 



^0(1,1,1) 



4i7r2 



ln=4, p2=_^2 



CI2 77 



TT 



(C.23) 



C.2.3 On-shell limit 



Now let us consider the on-shell limit pi = P2 = f^"^- For Ji(l,l,l), the two-fold 
parametric integral ( C.13 ) yields 



•/i(l,l,l)lp-,2_. = -ivr"/2r(3-f) / /- 

^1 



d^ d?7 



(1 + ^ + riY'^ [m?ri'^ - ^pl] 
Using Mellin-Barnes contour integral for the second denominator, we find 



2l3-n./2' 



(C.24) 



-i7^"/2f^2^n/2-3 



1 



X 



r(n-3) 27ri 

— ICXJ 

r(-s) r(n - 5 - 2s) r2(i + s) r (3 - f + s) . (c.25) 
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Closing the contour to the right, we get the result in terms of 2-^1 functions of the 
argument p|/(4m^). Note that the integral Ji(l,l,l) is not divergent in the on-shell 
limit, so that we can put n = 4 in Eq. ( C.24|) . For instance, the following simple 
representation in terms of CI2 can be mentioned 



"/l(l,l,l)ln= 



2m^ 



:4, Pj=p|=m2 



[Am? 



{CI 



2 ITT 



(C.26) 



with 



V'os = arctan J {Am? - pi) /pi, 



< P3 < 4m^ . 



Let us also consider the integral ^2(1; 1? 1) in the on-shell limit. Starting from the 
representation ( |C.14| ) in arbitrary dimension n and integrating over 1^, we get 



-^2(1, 1, I)|p2=p2^^2 



iivr"/^r (2 - 1) / - 



d?7 



2i3-n/2 



(C.27) 



It is easy to show that this three-point function (in the on-shell limit) reduces to a 
two-point function (|C.3|) with a shifted space-time dimension n — > n — 2, 



TT 



n 



-^2(1,1,0) 



n— >n— 2 



(C.28) 



Using Eq. 
^2(1,1,1 



we get 



ii^n/2 (^2)n/2-3r('2 



1, 3 -72/2 

3/2 



4^2 



(C.29) 



In the limit n ^ 4, because of the singular factor in front of the r.h.s. of Eq. (|C.28|) , we 
need to expand the two-point function up to the e term. This is how dilogarithms (or 
Clausen functions) arise in the finite part of the three-point function. 

The same result (|C.29| ) can be obtained from Eqs. (32) and (34) of |2^, taking into 
account that the first two arguments of the $2 function are Zi = Z2 = z = 1. Therefore, 
the sum over j in Eq. (34) of (corresponding to an 2-F1 function of unit argument) 
can be performed in terms of F functions. As a result, we arrive at Eq. ( |C.29|) . 

Moreover, using Kummer relations for contiguous 2-^1 functions, we get 



(n-3) 2F1 
Therefore, 



1, 2-n/2 
3/2 



{n-A){l-z) 2F1 



1, 3 -n/2 
3/2 



z\+l. 



(C.30) 



</2(l, 1, I)|p2^p2^^2 



2(n-3) 



{n — 4)(4m2 — p|) 



J2(l,l,0)- J2(0,l,l)b=„2 



(C.31) 
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Appendix D: Results for the x functions 

The one-loop results for Xi functions from Eq. ( p.25|) are collected below: 



X2\pI,pI,pI) 



X3\p1.p1.pI) 



[2 + {n- - m')^^ - (n - 2)ip{pipi^^ 



+ {n - 2)S,Mi[{piP2) + m]ipi- p^[n{pip2) + 2m jy^i 

+ {n - 3)^ m^pg^cs + p\{pI - m^)tii,i + m^{p\ - pi- m^)/ti,2 

-(2 - i)Mi{KQ^^ + «:i,2) + iMiK^,i + {n- 2)iim?\{p^p2) + im?]K^ , (D.l) 
g'^rj m Ca 



(47r)"/2 SICMi l^^ ~ OMi{p2P3)[iPiP2) + m]ipi- Kp^^i 
+ [2 + {n- 3)e]/CA^iV9i + - ^)^pI\{piP2) + m^] [p^^s ^ ^2^^^^^^ 

+ (2 - [(PlP2)/tl,l +P2'«1,2 + (P2P3)'«0,3 

+ - ?>)^pI [pIpI + m'^{pip2)\ /to,3 + (n - 3)^pl [pUpiPs] + m'^{,P2Pz) 
+ {n- ?>)ipl p\{p2Pz) + m^iPiPs) 1^1,2 + {n- 2)^m^/C R 



(D.2) 



^'"^ 1 [2 + (n - 3)e]/CA<i^i - (2 - 0-Mi(PiP3)[(PiP2) + m^]^! 



+e [(n - 4)/C - (n - 3)A^i] _ - 3)^plpl[ip,p2) + m^Ws 



(n - 3)^pI pI{piP2) + {pm){pl - rv? 



1^1,1 



-{n - 3)ipl{p2P'i)[{piP2) + m^]i^i^2 + ^P?/Cfi;i,2 

-(2 - i)Mi [Pi^l,! + (PlP2)/tl,2 + (PlP3)/«0,3 

+ (n - 2)em2/C(p^)-i[2(pip2) + m\K^^2 - 



(D.3) 



4)^m'lCpiipi - m'Mi[2{piP3) + ^(^2^3)]^! 



[2 + (n - 3)C]A^ip?(p2P3)¥'i -{n- 2)iMi{piP2){p2P3)Vi 
Ml[2{piP3)+i{p2P3)]K0,3 - Mi[2pl + i{pip2)]Ki^i - Ml[2{piP2)+ipl]Kl,2 



+ {n - 3)ipl [pIpI + m'^{pip2)\ Ko,3 + {n- 3)^pl [pUpm) + ^^(^2^3) 
+ {n - 3)^pI pj{p2P3) + m^{pip-i) Ki,2 + {n- 2)^m^/C k 



(D.4) 
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Appendix E: Transverse functions in an arbitrary gauge 

We here collect results for the transverse parts of the vertex, Tj, valid for arbitrary 
covariant gauge and dimension. 



E.l Transverse functions of diagram a 

Using the decomposition ( p.l5|) , all r's of diagram a can be presented as0 

A^"'\p1^p1^pI) = ^ ^ ^(4^)n/2 {4o"V2 + 4y [iPm)ni,l + iP2P3)ni,2+plK2,: 

+ti]2^ + 1^1,2 - 2^2,3) + tff + /tl,2 - 2k) 

+t^::H-.,i+-.,2)+&^^]. (E.l) 

^ Pi -pi > 

The results for the scalar functions t'^""* (which depend on the invariants ^1,^2)^3) ^-^^ 
listed below, for all eight r's: 

4y = -(^ - i)^^'\P\P2) - "ir?\ , 

f^l,! — '^1,3 — ''1,4 — U ) 

t'iS = \{n~i)mlC-\p,-p2)\ 

4y = -i/C-i{2(n - A)^M-,YM - "^')(P2 - rn')[(piP2) - m"] + (1 + OP3 " 

+ {n - l)lC-^pl[{piP2) - m\{l + OiPm) - (1 - 0"^'] + 4e(piP2)}, 
4;"^ = -i(^ - 1)^"'[(1 + 0(PiP2) - (1 - 0^1 , 

= ie/C-M(^ - 2) - (n - ^WM^\pI - plf] , 
4^3^ = -|(^ - 2)m2/C-H(l - i){plpl)-\piP2) + 2eA^2 V3[(PiP2) - m^] + AiKM-,^} , 
4!? = -|(^ - 4)/C-iA^2 H(l - 0^2 + 2ipl\plpl - m\p,p2)] - Aim'lC} , 
4^5^ = -\^'\p\plr'{{n - 2)(1 - 0(P1 - P2)'[P?P^ - m\p,p2)] + 2(n - 2)(1 - Om^/C 
+2(n - 3)em2A^2 - p2)'(p? - plf} , 



^^As mentioned in Section 3.2, the combinations of k's in Eqs. (E.l) and ( |E.2| ) are linearly dependent. 
As a result, we can (simultaneously) shift i^.i ti,i + 2ci, ti_2 ij.2 +P3Ci and t^^s — > ti^^ + {p\ —p\Yci 
(where Ci are arbitrary functions of the momenta, which can be chosen separately for each r^), without 
change of the value of the corresponding t;. This possibility allows us to write some of the t's in a more 
compact form. 
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IIC-'UpI+pI - 2m^)[(l + 0{pI+pI) - 2(1 - CW] - 4(rz - 2)/C 

+ (n - 1)}C-\pI - pl)'[{piP2) - m^] [(1 + 0(PiP2) - (1 - 0"i'] 
-2(n - 4)em^A^2 '(P? -P2)'[(PiP2) - m']{pl+pl - 2m^)}, 
0, 

3^/C-2{2(n - ^)im^lCM,\p\-plf{pl+pl - 2m2) 
-[(n - IM - Vlf + 4/C][(l + 0(PiP2) - (1 - 0^']} , 

^(n - 2)m2/C-H(l - 0»2)-M(Pi - P2)'(PiP2) - 2/C] + 2^M:,\pl - plf[{pm) 

^(n - 4)/C-H(l + - P2)' + 2im^M,\pl - plf[{pm) - m']} , 

^/C-2(p2 _ p2)2|^^ _ ;l ^ 20 [2/C - (piP2)(pi - P2f] - (n - < + 6e)/C 

-(1 - Om\n - 2)}C(pIpI)-\p,P2) - (n - l)(p, - p^)^ 

+{n - 3)^M^\p, - P2f[2m'plplpl - pUpmMpl + m') - 4m^/C]} , 

-l^mlC-^{{n - l)lC~V^[{p^p2) - m^] - 2{n - 5) 

+2(n - 4)}CM,\pI + pI- 2^2) - (n - 4)M,\pl - plfUpm) - m^]} , 

iem;C-^[(n - 3)M^\pI+pI - 2m') - {n - 1)IC-'] , 

0, 

i(n - 2)^m^]C-^M^' {pIpI)-'[2}C{pI+pI) - pIpIpI + ni?pl{p^P2)] , 
-\{n - 4)^m}C-'M^V3[{PiP2) - m'] , 

-\^m{(n - 2)K:-'[{plpl)-\pm){pl+pl) - 2] + 4(n - S)M^\pl+pl - 2m')} , 

->/C-H4(n - 4 - 0/C + (n - A)^}CM^\pl - Plfipl - ^rn') 
-(n - 3)e(p? - pD' - 2^pI[{p,P2) - m^]} , 
\imlC-\{n - 3)M^\pi - pD'iipm) - m'] + 2} , 
0, 

|(n - 2)em3A^2 i(p2p2)-i[2(p2 _p2)2^^^^^) ^ ^pjplpl - m'plipj+pl)] , 
|(n - 4)emAl2 Vsb? +P2 - 2m2) , 

l^m{pl - pl)'{{n - 2){pIpI)-' - 4(n - S)^!^ M(PiP2) - m']} , 
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C = -1}C-'M,\pI - pI){M2[{1 + OpI - 4^2] + 6^my,[{p,p2) - m']{pl + pi - 2m') 

+ in - l)/C-VI[(PiP2) - mT[-M2 + ipl{p\p\ - m%} , 
4;i^ = \Kr\p\-p\){2{n - 3)^m'M,\pl+pl - 2m') - {n - 1)1C-\{1 + i){p,p,) - (1 - 0^^]} , 

4;? = , 

4;3^ = -3^(n-2)m2/C-i(j9?-p^){(l-e)(p?P^)-^(piP2) + 2e^2Vi[(PiP2)-m2]} , 

4? = -U^-^)K,-'M~,\pl-pl){{l-i)M2 + 2ip\plpl-2im'[AK+pl{p,p,^^^ , 

45^ = j-,lC-\pl - pl){{n - 2)(1 - i)m'[{plpl)-\p,p,){pl+pl) - 2] - (n - 2)(1 + - P2)' 

4? = _ 4)A^2 + pI- 2m') - {n - 3)/C-i] , 

4;f^ = -|(^ - 3)em^2 H/C-V^[(piP2) - m'] + 2} , 
45^ = , 

4;3^ = -\{n- 2)im'M-,\plpl)-'{pl[2{p,p,) - m'\ + 4/C} , 

4;? = -i(^-4)e"^-^2-VL 

4? = iem{4(n - 3)^2 i[2/C + j9^(piP2) - m'p'^ - (n - 2)(p?p2)^i(p2 ^ ^2)^ ^ 

4? = - 6 - (n - 4)e]/C-V3[(PiP2) - m^] , 
4? = -i[^-6-(n-4)e]/C-S 

'-8,2 — ^^8,3 — '-8,4 — '^8,5 — U • 



E.2 Transverse functions of diagram h 

By analogy with diagram a, ( [E.l|) , all r's of diagram h can be presented as 



'^^^''^(P?,P2.P3) = \^^\^2 {4,0'Vl + 4,? [(PlP3)ft^l.l + (P2P3)/tl,2 +P3'«0,3] 



+4,2 (/^l,! + ^^1,2 - 2^0,3) + 4,3 + '^1,2 - 2k) 




(E.2) 
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The results for the scalar functions t\ ' (which depend on the invariants pi, pi, pi) are 
hsted below, for all eight r's: 

4y = lm}C-'M^^{{n - 4)(n - 6)^^m^lCM^'pl[{pm) + m^] + {n - A)ilCpl 
+ (n - 4)e'[(piP2) + m\{n - 2)K: + (n - 3)pl{p^p2)] 
-{n - 3)e(l - OM,pI - 2{2n - 2 - OMiUpm) + m^]} , 
4? = -|(^ - 3)em/C-^A4r'{[l -{n- SMpm) + [l-^-{n- 6)^/CMr>'} , 
4? = i^m[2(2n - 2 - - (n - 2)(n - 3)f A^r'] , 

4!? = ^(n - 2)^m{{plpl)-' - MT' -(n- 6)^m'M-['[(piP2) + m']} , 
tiA = - m^M-,'{?,i - 1 + (n - ^M-ApIpI + m\p,p^)]} , 
= ^m/C-^[2(2n - 2 - 0(Pi - {n - 2)^}C{plpl)-' {pi + pi)] , 

C = i^^"H[2 - e - (n - 3)e>^ - (n - 3)(n - 4)e^m2A1rV^[(PiP2) + m'] 

+(n - 1)(2 - 0/C-V^[(piP2) + mr -{n-l){n- S)e}C-Vs{PiP2)[(piP2) + m'] 
+[4(n - 3) - 2(2n - 5)^ - (n - 2)(n - 3)f ][(piP2) + m^]}, 

4? = Te^-'Hn - 1)(2 - 0[(PiP2) + m'] - {n - 3)e^[(n - l){p,p,) + (n - 3)m2/CA<^^]}, 

4? = -^/C-^[4(n - 3) - 2(2n - 5)^ -{n-2){n- 3)^^ , 

= - 2)w'K-\{n - 3)eM^'[{pm) + m'] - (2 - Oiplpl)-\PiP2)}, 

4? - 4)/C-^A1rM(2 - OA^i - - 3)e'biP2 + m2(piP2)]} , 

4!? = i/C-H(n - 2)(2 - OmWiPl)-\pm)ipl+pl) - 2] 
-[2(n-4)-(3n-8)e](pi-p2)'}, 

4y = ^/C-H(n-4)em^AtrM(PiP2)+m^][(n-3)e(p?-p^)^-4(4-0/C] 

+e[4 + (n - 4)e][(piP2) + m^][{n - 3){p, - P2)' - 4(n - 4)}Cipl)-'] 

+ (2 - 0[{n - 2){pI - pI)'{1 - IC-'M,) - IC-'Miplipi - P2)'] + 8(n-3)^/C 

+{n-3)e{pl-pl)' + {n-3)e{PiP2l{piP2)+m'][^ + {n-l)]C-\pl-pl)']}, 

4? = ^,IC-\pl)-'{(n-3)Cm'M^'[(n-3)apl-pl)'-^(^-m 

+ [{n - 1)}C-\pI - plf + 4][(n - 3)e{piP2) - (2 - 0(PiP2) - (2 - 0^']}, 
4;2^ = -^e/C-^(pD-M4+(n-4)e][4/C+(n-3)(p?-p^)2] , 
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^(n - 2)m'lC-\pl)-'{CMT'[{piP2) + m^][4(4 - 0/C -{n- S^pj - plf] 

+(2 - OiplplrVslipmKPi - P2? - 2/c]} , 

^{n- A)lC-\pl)-'{{n - 3)eM^\pl - plf[plpl + m\pm)] 
-(2 - i)plipi - P2f - m - 0^-Mr'biP2 + ^'(PiP2)]} , 

j-,lC-\pl)-\2 Oipl - Plf{<n - 2WlC{p\pl)-' ~ 2{n ~ 3)^{p, - p,f 
-{n - 2){p^-p2f[m\plpl)-\pm) + 1]} , 

\imlC-^M^^{{n - 4)(n - Q)im'^lCM^^pl + (n - 2)(n - 4)^/0 

+ [1 + (n - ^)i]pl{piP2)[n - 4 - (n - l)lC~^Mi] - {n - l)m^K,-^Mipl 

+ (n - 4)[1 - (n - 3)^]m2p2 _ _ 2)(n - 3)CMi}, 

-(n - l);C-^A<i + (n - 3)(1 + 3^)} , 

^(n - 2)(n - 3)^2^/C-iMr'[(PiP2) + m^] , 

^(n - 2)im1C-^M^\Mi{piP2){p\pl)-^ - (P1P2) - + (n - 3)em' - (n - 6)em2/CAt 
-3^(n - A)^mlC-^M^^{{piP2) +m^ + {n- 3)e(piP2) + (n - 6)em'/CA<r'} , 
-^(n - 2)em/C-H(p?P^)-^(piP2)(p? - 2] , 

-\mK:-'{ipl\pl+pl + (2 - e)m2][(n - A)1CM-,' -{n- 3)] 

+ [(n - 3)(2 - - lM[<^iP2) + m^] + [12 - 2^ + (n - 4)^(4 - 0]/C}, 

j^(n-3)em/C-^A^rH2/C(p?+p2_2rri2) _ (p2-p2)2[^^^^^)^^2] _ ^^2^2^2^^2^^^^^)|^^ 

^(n - 2)emA^rH(2 - - [MMpIY' - +P2)} , 

3^(n - 4)emA1^i[p2 +p2 + (2 - ^m^] , 

j-^^mK:-'M-,\pl - plf{{n - 3)(2 - ^Mpl + /^^(pip^)] 
+(n - 2)1CMMpI)-' -{n- 3)(p? +p^)[(piP2) + m^]} , 

-1,1C-\pI - pl){{n -3)(n- A)e(m'MT'pl + l)[(piP2) + m'] 

+in - l)(n - 3)e}C-'plipiP2)[iPiP2) + m'] + [{n - 2)(2 - + (n - 3)e]pl 
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-2[2 - (n - SmiPm) + m'] - (n - 1)(2 - O^C-'Mipl}, 

C = i/C-2(p^P^){(n-l)(2-0[(piP2)+mV(n-l)(n-3)C^(piP2)-(n-3)^eWK:A<^i}, 
4!? = -^^-'(P? -pD[4 - 2{n - 3)e - (n - 3)(n - 4)^^], 

C = ^(^ - m-\pi - pDiin - 3)eWA1rM(PiP2) + m'] - (2 - 0^'(P?P2)-'(PiP2)}, 
C = U^-^)IC-'M^\pi-pl){{2-0M,-in-3)e\plpl + m'{p^p2)]}, 

C = ^/C-^(p? - 2)(2 - 0mWiPlr\PiP2){pl+pl) - 2] 

+ - P2)'[(n - 4)(2 + - 2(n - 3)^']} , 

4? = - 4)Mr' - - 3)/c-^] , 

^(16) _ n 

4;? = 3^(n-2)em[(p?p^)-i-A1r^ 

4? = -]^(^-4)emAir', 

4!? = -i^(^-2)Mp?p^)-^(p?+p^), 

4? = |/C-H(n - 4)e(4 - OiCirn'MTVs + l) + {n- 3)^(4 - 0p^(PiP2) 
+(6-0[P3(PiP2)+m2p^ + 2/C]}, 

4? = |/C-i>trH(^ - 3)e(4 - 0b?P2 + m'ipm)] + (6 - OA^i) , 

.(16) _ ,(16) _ p, 

4;? = -i^(n-2)e(4-eK-MrS 
4y--i^(^-4)e(4-eKA^r'- 

In the hmit m — > 0, r^^^ r]^^ Tg^-*, and Ty^^ vanish. We also note that all 4,4^^ and 

4,4^'* are proportional to (n — 4), as they should, since the transverse part cannot contain 
UV-poles in e at one loop. 
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